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Figure 

1 Model  of  a cable  non-unifonnity  as  a symmetrical  T section. 

The  network  is  to  be  loaded  on  the  line  at  the  location  of 
the  non-uniformity.  The  lumped  impedance  and  admittance 

are  to  be  calculated  from  two  quasi-static  boundary-value 
problems. 

2 Section  of  a uniform  transmission  line  loaded  at  discrete 
points  with  four-terminal  networks. 

3 Graphical  solution  of  the  dispersion  relation  (42)  of  a 
periodically  loaded  transmission  line  in  the  (cos  ka,  m) 
plane.  The  frequency  ranges  for  which  -1  i cos  ka  s l 
correspond  to  the  passbands  (dotted  areas). 

4 Plot  of  frequency  m as  a function  of  propagation  constant  22 

k , showing  the  passband-stopband  structures.  The  curves 

£or  negative  k and  u can  be  obtained  by  reflections 
with  respect  to  the  horizontal  and  vertical  coordinate 
axes. 
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SECTION  I 


INTRODUCTION 


At  sufficiently  low  frequencies  when  only  the  TEM  mode  is  dominant,  the 
propagation  of  voltage  and  current  disturbances  along  a uniform  cable  can  be 
described  by  the  pair  of  transmission  line  equations: 


3x 


(1) 


31 

3x 


(2) 


For  a time-harmonic  excitation  with  time  factor  exp (jut)  , these  equations 
become 


dx 


ZI 


(3) 


ii 

dx 


- YV 


(4) 


where 


Z = juL  + R 

(5) 

Y = juC  + G 

(6) 

are  the  impedance  and  admittance  per  unit  length  of  the  uniform  line  at 
angular  frequency  u • Excitations  with  general  time  variation  can  be  resolved 
into  time-harmonic  components  by  the  method  of  integral  transform.  Consequently 
their  treatment  can  also  be  based  on  equations  (3)  and  (4). 

The  uniform  line  is  of  course  an  Idealization.  In  practice,  non-uniformities 
occur  in  the  shape  of  cable  bends,  shield  defects,  or  nearby  conductors.  These 
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produce  local  deviations  of  Y and  Z from  their  uniform  values,  resulting 
in  a scattering  of  the  waves  in  the  cable.  If  the  dimensions  of  the  non- 
uniformities are  small  compared  to  a wavelength,  the  scattering  effects 
can  be  described  by  an  equivalent  point  loading  of  the  cable.  Fig.  1 
shows  a model  of  the  non-uniformity  as  a symmetrical  T section.  The 
network  is  assumed  to  have  zero  spatial  extension,  and  is  to  be  inserted 
into  the  uniform  line  at  the  location  of  the  non-uniformity.  The  lumped 
impedance  and  admittance  elements  and  Y^  of  the  discontinuity  are 

to  be  calculated  from  appropriate  quasi-static  boundary-value  problems. 

The  general  situation  of  a transmission  line  of  finite  length,  driven  at 
one  end  by  a harmonic  voltage  or  current  source,  terminated  at  the  other 
by  a given  load,  and  loaded  at  an  arbitrary  configuration  of  points  in 
between,  can  be  analyzed  by  the  powerful  transmission  matrix  method  to  be 
described  below  in  Section  III. 

Of  special  interest  to  aircraft  EMP  internal  coupling  studies  is  the 
analysis  of  wave  propagation  along  a periodically  loaded  transmission  line. 

Here  the  loads  are  Identical  in  structure,  and  are  inserted  along  the  line 
at  regular  intervals.  Examples  of  such  periodic  loading  in  aircraft  cables 
can  be  found  in  periodically  applied  cable  clamps  and  periodic  airframe 
members  over  which  the  cable  runs  are  anchored.  It  is  well  known  that  a 
periodically  loaded  line  acts  like  a bandpass  filter.  The  line  transmits 
waves  with  frequencies  lying  within  certain  discrete  bands,  while  waves 
with  other  frequencies  are  effectively  stopped.  It  can  therefore  restrict 
the  range  of  EMP  energy  that  can  be  delivered  through  it  to  the  load.  Hence 
a determination  of  the  passband-stopband  structures  as  a function  of  the 
loading  is  extremely  helpful  to  the  assessment  of  EMP  vulnerability  of  air- 
craft systems. 

The  objective  of  this  report  is  to  calculate  the  propagation  characteristics 
of  an  infinitely  long  transmission  line,  loaded  periodically  with  four-terminal 
networks  of  the  type  shown  in  Fig.  1.  This  study  will  consist  of  deriving 
the  so-called  dispersion  relation  connecting  the  frequency  with  the  propagation 
constant,  showing  its  explicit  dependence  on  the  loading  period,  and  on  the 
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distributed  and  lumped  impedances  and  admittances  of  the  uniform  line  and  its 
lumped  loads.  The  dispersion  relation  governs  the  propagation  of  the  eigenmodes 
of  the  infinite  line.  The  reason  for  the  interest  in  these  infinite-line  modes 
is  that  the  voltage  and  current  excitations  on  the  finite  line  encountered  in 
practice  are  composed  of  these  eigenmodes.  The  composition  is  determined  by 
taking  the  applied  source  and  the  load  termination  as  boundary  conditions  at 
the  two  ends. 

The  contents  of  this  report  are  distributed  as  follows.  Section  II 
summarizes  the  salient  features  of  waves  on  a periodic  line.  The  transmission 
matrix  method  for  calculating  wave  propagation  along  transmission  lines  is 
introduced  in  Section  III.  In  Sections  IV,  V and  VI  the  dispersion  relation, 
the  passband-stopband  structures  and  the  phase  and  group  velocities  of  a 
transmission  line  loaded  periodically  with  symmetrical  T sections  are 
explicitly  derived.  The  report  concludes  with  an  illustrative  example 
worked  out  in  detail  in  Section  VII. 


SECTION  II 


GENERAL  PROPERTIES  OF  WAVES  ON  A PERIODIC  LINE 


The  phenomenon  of  wave  propagation  in  periodic  structures  occurs  in  many 
areas  of  physics  and  engineering.  Over  the  years  its  principal  features 
have  been  delineated  and  the  general  procedures  for  its  study  worked  out. 
Before  embarking  on  the  specific  calculations  for  the  periodically  loaded 
transmission  line,  it  is  advantageous  to  anticipate  the  results  by  summarizing 
below  the  general  properties  of  waves  in  a one-dimensional  periodic  structure. 

1.  At  a given  frequency  u , the  waves  on  a line  of  period  a have 
Che  general  form 


■Hx)  = e 


(7) 


where  k and  u,.  (x)  are  determined  bv  the  detailed  structure  of  the  line,  and 
Uj^(x)  is  a periodic  function  with  period  a . Thus  Che  excitations  are  "plane 
waves"  modulated  by  the  periodicity  of  the  line.  This  result  is  known  in 
mathematics  as  Floquet's  theorem. 


2.  Upon  analyzing  the  periodic  function  into  a Fourier  series 


u,(x) 


exp 


jlnTx^ 


(8) 


one  concludes  from  (7)  Chat  the  excitations  consist  of  an  infinite  sequence  of 
spatial  harmonics  with  propagation  constants 

kn  = k + -^  , n = 0,  ±1,  ±2,  ...  (9) 

Hence  at  each  frequency,  there  are  an  infinite  number  of  eigenmodes. 
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In  fact,  the 


3.  The  frequency  oi  is  an  even  periodic  function  of  k . 
dispersion  relation  is  of  the  general  form 

cos  ka  = F(u)  (10) 

4.  When  there  are  no  resistive  elements  in  the  line,  there  exist  an 
infinite  series  of  frequency  bands  in  which  k is  alternately  real  and  purely 
imaginary.  The  line  is  a bandpass  filter. 

These  general  results  will  be  verified  and  illustrated  by  the  calculations 
in  Che  following  sections. 


SECTION  III 


"I 


TRANSMISSION  MATRICES 

Consider  a length  of  uniform  transmission  line  containing  a number  of 
point  loads,  as  shown  in  Fig.  2.  Each  load  is  assumed  to  have  the  structure 
shown  in  Fig.  1.  The  variation  of  the  voltage  and  current  along  the  uniform 
sections  is  described  by  the  transmission  line  equations  (3)  and  (4).  The 
variation  in  the  loads  is  governed  by  Kirchhoff's  laws. 

At  any  reference  plane  along  the  line,  such  as  plane  A in  Fig.  2,  the 
voltage  and  current  can  be  combined  into  a two-component  vector: 


Then  the  voltage  and  current  at  any  other  reference  plane  in  the  same  uniform 
section  of  the  line,  such  as  plane  B in  Fig.  2,  are  related  to  those  at 
plane  A through  a transmission  matrix  T^: 


= T 


^V(A)^ 


(11) 


This  matrix  is  determined  by  solving  equations  (3)  and  (4).  It  is  given  by 


'cosh  ya 


-Z  sinh  ya ' 
0 


u \ 1 . , , 

— sinh  ya  cosh  ya 

^0 


(12') 


where  a is  the  separation  between  planes  A and  B , and 


='/f  • ^ 


(13) 


are  respectively  the  characteristic  impedance  and  the  propagation  constant  of 


: v''  ■ 


sc 


r 


1 


the  uniform  section. 

Similarly,  the  variation  of  the  voltage  and  current  across  a four- 
terminal  network,  such  as  between  reference  planes  B and  C in  Fig.  2, 
can  be  described  by  a transmission  matrix  : 


V(C)\ 


1=  T 


,i(cV 


/V(B)  V 


(14) 


Upon  applying  Kirchhoff’s  laws  to  Fig.  1,  one  finds  that 


"d=‘ 


-Yd 


■^d  - 


1 + '5Y,Z 


d d 


(15) 


It  is  to  be  noted  that  both  T and  T,  have  determinant  1: 

u d 


det  T =1  , det  T,  = 1 

u d 


(16) 


Matrices  with  this  property  are  said  to  be  unimodular. 

The  relation  between  the  voltage  and  current  at  reference  plane  C and 
those  at  plane  A can  be  obtained  by  matrix  multiplication.  Thus 


V(C)\  /V(A)') 


1(C), 


j(A)y 


(17) 


wnere 


4T  T 
'11  ^12 


VT  T 

:i  2: 


T.T 
d u 
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(13) 


with 


i 


I 

r 


= (1  + JiY^Z^)cosh  Ya  + ^(1  + 5jY^Z^)sinh  ya 

0 

T««  = (1  + ^Y,Z,)cosh  ya  + Y.Z  sinh  ya 

u Q GO 

=■  -Z,(l  + itY.Z,)cosh  ya  - Z (1  + isY_,Z,)sinh  ya 
a d d 0 d d 

= -Y^cosh  ya  - -^(1  + 55Y^Z^)sinh  ya  (19) 

It  is  easy  to  show  that  T is  also  unimodular: 

det  T = Tj^^T22  - Tj^2'^21  " ^ 

Using  the  transmission  matrix  method,  the  calculation  of  voltage  and 
current  propagation  along  a finite  transmission  line  with  discrete  point  loads 
becomes  a simple  matter  of  2x2  matrix  multiplications. 


SECTION  IV 


DISPERSION  RELATION 


Suppose  in  Fig.  2 the  transmission  line  section  between  reference  planes 
A and  C is  repeated  an  infinite  number  of  times  to  the  right  and  to  the 
left.  One  obtains  in  this  way  an  infinite  periodically-loaded  transmission 
line  of  period  a . If,  at  a given  frequency  u , a wave  is  to  propagate 
down  this  structure,  the  voltage  and  current  in  reference  plane  C can 
differ  from  those  at  plane  A by  at  most  a phase  factor.  This  is  because 
planes  A and  C are  separated  by  exactly  one  period,  and  are  hence 
physically  equivalent.  Let  this  phase  factor  be  denoted  by  exp(-jka)  , 
that  is. 


(21) 


This  result  actually  follows  from  Floquet's  theorem  (7),  and  k can  be 
identified  with  the  propagation  constant  of  an  eigenmode.  Relation  (21) 
is  compatible  with  relation  (17)  only  if  the  phase  factor  is  an  eigenvalue 
of  the  transmission  matrix  T . One  therefore  requires  that 


0 


This  condition,  by  (20),  becomes 


(22) 


(Til  -22^®'^*^  +1-0 


(23) 


The  solution  is 

cos  ka  - -i(T^^  + T.,^) 


(24) 


Upon  substitution  of  the  matrix  elements  from  (19) , equation  (24)  becomes 


cos  ka  = (1  + iiY,Zj)cosh  ya  + + 

da  ZZ  d 


h'i.zl  + Y,zf)sinh  ya 


The  quantities  Y^  , Z^  , Z^  and  y are  all  functions  of  the  frequency  lo  . 
Therefore  (25)  is  a dispersion  relation  connecting  u with  the  propagation 
constant  k . It  agrees  with  the  general  form  (10) , If  k is  a solution 
of  (25)  at  a given  oi  , then  so  are  -k  and 


as  well  as  their  negatives  -k^  , because  of  the  evenness  and  periodicity  of 
the  cosine  function. 


SECTION  V 


PASSBANDS  AND  STOPBANDS 

In  many  practical  situations  the  resistive  elements  in  the  line  are 

negligible,  so  that  all  the  impedances  and  admittances  are  expressible  in 

terms  of  inductances  and  capacitances  L , C , L,  and  C,  : 

d d 

Y = juC  , Z = juL 

^d  = J“^d*  ^d  = J“^d  (27) 

From  these,  one  has 

Zq  =yi"  . Y = jw  (28) 

In  this  case,  relation  (25)  becomes 

2 , — 

cos  ka  = (1  - L ,C  ,)cos(o)a/LC) 
d d 

fF  2 

(1  - sin(wjai/L3)  (29) 

d 

The  right-hand  side  is  an  even  function  of  the  frequency.  Equation  (29)  can  be 
put  into  a simpler  form  by  introducing  an  amplitude  A(u)  and  a phase  >(u) 
such  that 


A(u)  cos  !(i(u)) 


1 - 


^d^d 


A(u)  sin  iji(u) 


(30) 


Then  (29)  becomes 


cos  ka  » A(u))cos(ua»'LC  + ij>(w)) 


(31) 


The  right-hand  side  is  a quasi-sinusoidal  function  with  a modulated  amplitude. 
After  some  algebra  it  is  found  that 


A^(ui) 


1 + 


•|l  - 


4L 


\ 


d^d 


(32) 


which  is  greater  than  1 . From  this  one  concludes  that  there  exist  an  infinite 
number  of  frequency  ranges  over  which  the  right-hand  side  of  (31)  is  between  the 
limits  1 and  -1  , and  other  ranges  over  which  the  expression  is  outside  those 
limits.  Equation  (31)  shows  that  k is  real  in  the  first  case,  implying 
propagation;  it  is  imaginary  in  the  second  case,  indicating  attenuation. 
Consequently  one  has  an  infinite  sequence  of  alternating  passbands  and  stop- 
bands.  The  band  boundaries  are  determined  by  the  condition 

cos  ka  » ±1  (33) 

or 

k = -^  , n = ±1,  ±2,  ...  (34) 

d 

This  condition  corresponds  to  a total  reflection  of  the  waves  at  the  loads. 
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SECTION  VI 


PHASE  VELOCITY  AND  GROUP  VELOCITY 

Let  the  right-hand  side  of  the  dispersion  relation  (29)  be  denoted  by 
F(cj)  : 


cos  ka  = F(ai)  (35) 

If,  at  a given  frequency  u in  a passband,  k is  a real  solution  of  (35), 
then  so  are  -k  and  the  members  of  the  infinite  sequence 


k + 


±1,  ±2, 


(36) 


as  well  as  their  negatives  -k^  . The  waves  on  the  transmission  line  therefore 
consist  of  an  infinite  number  of  spatial  harmonics  traveling  to  the  right  or  to 
the  left.  Each  harmonic  has  a phase  velocity 


V 


ph 


(37) 


which  is  different  for  different  n's  . The  ± sign  refers  to  the  two 
possible  directions  of  propagation.  Thus,  at  a given  frequency,  there  is  no 
unique  phase  velocity. 

There  is,  however,  a unique  group  velocity.  The  group  velocity  is  the 
velocity  at  which  energy  is  transmitted  by  the  waves,  and  is  a very  useful 
quantity  in  analyzing  propagation  characteristics.  It  is  given  by 


From  (35)  one  finds  that 


3u) 


f 2 

_ ^ a>l  - F (w) 
'^g  " 5F(u)) 

3(jj 


(38) 


(39) 
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which  has  a uniquely  defined  magnitude  at  a given  frequency.  Again  the  ± 
sign  refers  to  the  two  possible  directions  of  propagation.  By  (33)  and  (35) , 
F(a)  » ±1  at  the  band  boundaries.  Therefore  one  has  the  general  result  that 
the  group  velocity  vanishes  at  band  boundaries. 
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SECTION  VII 


ILLUSTRATIVE  EXAMPLE 


In  this  section  a numerical  example  will  be  worked  out  to  illustrate  the 
general  results  derived  above.  Consider  the  dispersion  relation  (29)  applicable 
to  the  case  of  negligible  resistive  circuit  elements.  It  is  convenient  to 
introduce  three  dimensionless  parameters 


= uai''LC 


a 


a^C 


3 


(41) 


Hence  0 serves  as  a dimensionless  frequency  variable,  and  a and  B are 
parameters  characterizing  the  periodic  loads.  Then  (29)  becomes 


cos  ka 


(42) 


Take  a specific  set  of  values  for  a and  B : 


01  - 0.2  , 3 =■  0.5  (43) 

This  choice  for  3 is  quite  typical  of  practical  situations.  The  choice  for 
a corresponds  to  a fairly  heavy  loading  of  Che  line.  It  will  be  seen  below 
that  the  loading  (43)  produces  at  low  frequencies  a reduction  of  the  group 
velocity  to  about  2/3  the  uniform  line  value. 

The  passband  and  stopband  structures  can  be  found  graphically  by  plotting 
out  the  right-hand  side  of  the  dispersion  relation  (42)  as  a function  of  the 
dimensionless  frequency  variable  iJ  , as  is  done  in  Fig.  3.  The  curve  oscillates 
with  ever  increasing  amplitude  as  fj  -►  ± ^ . Those  portions  of  the  curve  lying 


Figure  3.  Graphical  solution  of  the  dispersion  relation  (-12) 
of  a periodically  loaded  transmission  line  in  the 
(cos  ka,  u)  plane.  The  frequency  ranges  for  which 
-1  cos  ka  ^ 1 correspond  to  the  passbands 
(dotted  areas). 


between  the  horizontal  lines  cos  ka  = 1 and  -1  correspond  to  the  passbands. 
They  are  indicated  by  the  dotted  areas  in  the  figure. 

Each  intersection  of  the  curve  in  Fig.  3 with  a horizontal  line  yields  a 
pair  of  coordinates  (cos  ka  , u)  which  form  a solution  of  the  dispersion 
relation.  Each  value  of  cos  ka  corresponds  to  an  infinite  sequence  of 
propagation  constants  ±k^  according  to  (9).  In  a passband,  these  constants 
are  real.  To  study  the  phase  velocity  and  group  velocity,  it  is  convenient 
to  plot  the  solutions  in  the  (k  , u)  plane.  This  is  done  in  Fig.  4.  Each 
continuous  horizontal  curve  corresponds  to  a passband.  The  passbands  are 
separated  by  gaps  representing  the  stopbands.  It  is  seen  that,  as  the 
frequency  increases,  the  passbands  contract  while  the  stopbands  widen.  Thus 
most  of  the  low-frequency  waves  are  passed,  while  most  of  the  high-frequency 
ones  are  stopped.  This  is  not  surprising  since  four-terminal  networks  of 
the  type  shown  in  Fig.  1 function  as  a low-pass  filter  when  connected  in 
tandem. 

The  phase  velocity  at  frequency  m , as  defined  in  (37),  is  proportional 
to  the  slope  of  the  straight  line  drawn  from  the  origin  to  a point  with  ordi- 
nate (1)  on  a curve  of  Fig.  4.  On  the  other  hand,  the  group  velocity  is 
proportional  to  the  slope  of  the  curve  at  that  point.  From  Fig.  4,  it  is  seen 
that  the  group  velocity  is  greatest  at  low  frequencies.  One  can  find  the  low- 
frequency  limit  of  the  group  velocity  by  expanding  (42)  in  ascending  powers 
of  oj  . The  result  is 


k r nr  "1 

m = i ^ 1 + a (1  + S) 

/lc  L ^ J 
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Therefore,  by  (43),  the  group  velocity  at  low  frequencies  is 


3(0  _ ^ 0. 68 
“ 3k  ■ “ 


v'lC 


(45) 


showing  a reduction  of  some  32%  from  the  uniform  line  value  by  the  loading. 


Figure  4.  Plot  of  frequency  u as  a function  of  pro- 
pagation constant  k , showing  the  passband- 
stopband  structures.  The  curves  for  negative 
k and  uj  can  be  obtained  by  reflections  with 
respect  to  the  horizontal  and  vertical  coordinate 
axes . 


As  the  frequency  is  increased,  Che  passbands  become  progressively  narrower. 
In  Fig.  4,  they  are  seen  to  approach  the  limit  of  a horizontal  line.  The  group 
velocity  therefore  tends  to  zero,  implying  zero  energy  transmission.  In  this 
limit,  the  excitation  on  the  line  consists  of  standing  waves  in  the  uniform 
sections  of  the  line,  bouncing  back  and  forth  between  two  loads. 
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